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On the basis of Quantum Monte Carlo simulations of the two-dimensional Hubbard model which 
cover the doping range from the under- to the over-doped regime, we find that the single-particle 
spectral weight A(k, u) qualitatively reproduces both the momentum (d^.^-symmetry) and doping 
dependence of the pseudogap as found in photoemission experiments. The drastic doping dependence 
of the spin response Xs (?, w) which is sharp in both q (~ (tt,tt)) and u in the under-doped regime but 
broad and structureless otherwise, identifies remnants of the antiferromagnetic order as the driving 
mechanism behind the pseudogap and its evolution with doping. 

PACS numbers: 74.20.-z,74.72.-h,75.50.Ee,79.60.Bm 

Exciting progress in the microscopic understanding of 
high- Tc superconductors has recently come from the ob- 
servation of a normal-state pseudogap of order of the ex- 
change energy J [jj] and a lower energy excitation gap 
of the order of the superconducting gap Jl|-§|. Angle- 
resolved photoemission spectroscopy (ARPES) demon- 
strated that both high- and low-energy gaps are con- 
sistent with d x 2_ y 2 -symmetry Jj]-^|]. In addition, both 
gaps have a more or less identical doping dependence, 
which may be a key observation to unlocking the mys- 
tery of the cuprates: just below optimal hole concentra- 
tion, the centroids in the spectral weight A(k,ui) near 
(tt, 0) move to higher binding energy, and a portion of 
the large Fermi surface seems to disappear. These find- 
ings have been interpreted as the opening of a pseudogap 
with maximal energy J ~ 2007716^ near (tt, 0) [ID. Si- 
multaneously, a normal-state gap with energy ~ 20meV, 
inferred from the leading edge in A(k,u>), opens up in 
this under-doped regime. Both of these gaps vanish in 
the over-doped regime, and the superconducting gap also 
rapidly decreases |5| . This empirical correlation between 
the disappearance of the order J pseudogap and the de- 
crease of superconducting pairing strength suggests that 
the high-energy features at (tt, 0) are closely related to 
the pairing interaction Q. 

In this letter, we address the microscopic mechanism 
behind the opening of this pseudogap and its evolution 
from under-doped to over-doped regimes. We present 
Quantum Monte Carlo (QMC) simulation results on the 
two-dimensional Hubbard model with on-site interaction 
U = 8t which demonstrate that the single-particle spec- 
tral weight A(k, to) reproduces the most salient ARPES 
features as function of doping. In particular, the QMC 
data reproduce the momentum {d x 2_ y 2 -symmetry) and 
doping dependence of the pseudogap. 

Earlier finite-temperature QMC work on the Hub- 
bard model |j],H has produced results showing a quasi- 
particle-like band with a dispersion below the Fermi level 
that is essentially unaffected by doping. On the other 



hand, groundstate exact diagonalizations [p[]lOf| of the 
two-dimensional t — J model for small clusters around 
optimal doping find a signal in the spectral weight cor- 
responding to an insulator-like "shadow" structure; at 
larger doping this signal vanishes. We present in this 
work first data of the spectral weight A(k, ui) obtained by 
Maximum-Entropy techniques for previously inaccessible 
temperatures (T = 0.25t), showing that if the tempera- 
ture in the dynamical QMC simulation is lowered below 
a threshold temperature T* ~ 0.3i in the under-doped 
regime the quasi-particle band is substantially deformed, 
resulting in the opening of a pseudogap. 

Our results provide numerical evidence of earlier the- 
oretical conjectures that the deformation and the pseu- 
dogap are intimately related to antiferromagnetic spin 
fluctuations JTl]— p~^] . This is demonstrated by studies of 
the two-particle excitation spectrum over a wide range 
of dopings from under-doped to over-doped. They reveal 
that it is the strongly doping-dependent spin response 
and not the essentially doping-independent charge re- 
sponse which tracks the doping dependence of the pseu- 
dogap. The spin response, \s (<f, w), displays a sharp 
structure in Im Xs at wave-vector q = Q = (tt, tt) and 
at energy uj(Q) = lu*. Not only near half-filling but also 
up to optimal doping the spectral weight is distributed 
around an energy dispersion which still closely follows the 



spin- wave dispersion w(q) 



,SDW 



(q), where oj sdw is 



calculated within the spin-density-wave approximation. 
For temperatures T <T* , the antiferromagnetic correla- 
tion length £(T) becomes larger than the lattice spacing 
ao and, as a consequence, the quasi-particle is strongly 
dressed by spin fluctuations. As soon as the system en- 
ters the over-doped regime, the spin response is no longer 
sharply peaked near Q = (tt, tt) and uj* . It spreads in en- 
ergy by an order of magnitude (the scale changes from J 
to Ekin ~ 8£) , and is accompanied by a similar change of 
the scale of the bandwidth for single-particle excitations. 
Thus it is this unique doping dependence of the magnetic 
excitations which establish the pseudogap as being due 
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to antiferromagnetic spin fluctuations. 

The single-band two-dimensional Hubbard model has 
the standard Hamiltonian 

-pi^(n lT +n a ) (1) 

i 

on a square lattice, where t is the nearest-neighbor hop- 
ping, Cj i<T destroys a particle of spin a on site i, and 
n ia = c irr Ci t(T . The chemical potential u sets the fill- 
ing (n) — (nif + riii). The spectral weight A(k,uj) 
is inferred from high-quality QMC data by apply- 
ing state-of-the-art "Maximum-Entropy" techniques [ 15) . 
This method has previously been used to resolve a quasi- 
particle-like dispersive band of width J at half-filling and 
to follow its evolution from the insulator to the metal [0] . 

In Fig. 0(a) we compare the QMC results for the 
single-particle spectral weight A(k,u>) for the under- 
doped (1% and 5%) and maximally doped (13%) regimes 
with the ARPES data from Ref. @ in Fig. 0(b). The 
theoretical results were obtained for temperatures of 
T = 0.33t and T = 0.25t, on-site Coulomb repulsion 
U = 8t and 8x8 lattices. To facilitate detailed compari- 
son with Fig. 0(b) of D. S. Marshall, Z.-X. Shen et al. §, 
we use in Fig. 0(a) the same u vs k "band structure" plot 
with all spectra referring to the same chemical potential 
fi(u = 0). 

In accordance with the ARPES data, we find a dra- 
matic change in the electronic structure near (n, 0) with 
doping. In the under-doped (1% and 5%) regime, the 
features near (tt, 0) move to lower binding energy as dop- 
ing is increased. At about 13% doping || (full circle in 
Fig. 0(a)) the pseudogap vanishes. 

In addition, a portion of the large Fermi surface, which 
is closed around (tt, tt) at T — 0.33i (i.e. has Fermi 
level crossings between (tt, 0) and (tt, tt) as well as be- 
tween (tt, tt) and (0,0)), seems to disappear at the lower 
temperature, T = 0.25i. This is indicated in Fig. 0(a) 
by the downturn of the quasi-particlc-like band between 
(tt, 0) and (tt, tt). In the experiment, this behavior in the 
under-doped regime has been interpreted as the opening 
of a pseudogap in the underlying Fermi surface near the 
(tt, 0) to (tt, tt) line [Q. The momentum dependence of the 
pseudogap is also in agreement with experiment in that 
its largest value, which is of the order of the exchange 
coupling J ~ At 2 j\] , occurs near (7r,0) and in that it 
essentially vanishes along the (0, 0) to (tt, tt) line. This 
behavior is consistent with d 2 .2_ y 2-symmetry. At 13% 
doping we obtain a Fermi surface crossing near (tt, 0) 
(full circle in Fig. 0(a)) independent of the two chosen 
temperatures. Thus, at optimal hole doping the Fermi 
surface is large, consistent with both ARPES 0-0] and 
earlier QMC calculations H. 
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FIG. 1. The dispersion of the peaks in the single-particle 
spectral weight from (a) QMC simulations of the Hubbard 
model at densities ranging from the under-doped to the maxi- 
mally doped regime where peaks in A(k, u) axe represented by 
error bars and (b) peak centroids in A(k, ui) from under-doped 
and maximally doped ARPES experiments after Ref. [1]. 

In the QMC simulation, the opening of a pseudogap 
in the under-doped regime shows up below a "crossover" 
temperature T* ~ 0.3£. To illustrate this crossover, we 
consider in Fig. 0(a) - 0(b) the role played by temper- 
ature and its influence on the spectral weight distribu- 
tion in more detail. We observe two general features 
in both spectra: a (several t) broad "incoherent back- 
ground" both below and above the Fermi surface and a 
dispersing structure with a smaller width of order of a 
few J around u — u. In our earlier work 0] in which 
the lowest temperature accessible was T = 0.33t, this 
dispersing quasi-particle-like band was shown to have its 
maximum at (tt, tt). In this case, which is reproduced 
in Fig. 0(a), there is a "large" Fermi surface centered 
around (tt, tt). However, by lowering the temperature to 
T = 0.2bt, the structure forming the maximum at (tt, tt) 
looses weight and the new valence-band maximum seems 
to be shifted to (tt/2,tt/2) or (tt, 0) (Fig. |(b)) |§. In 
agreement with the ARPES data in Fig. 0(b), we observe 
the downturn ("shadow" structure) as well as a drastic 
(by about a factor of 10) spectral weight loss when fol- 
lowing the quasi-particle band from (tt, 0) to (tt, tt). It is 
well known [0 that inclusion of higher-neighbor (t', etc.) 
interactions lifts the degeneracy of the points (tt/2,tt/2) 
and (tt,0) at half-filling and pushes the (tt/2,tt/2) point 
up in energy. This opens up the possibility of a small 
Fermi surface ("hole pockets") around (tt/2,tt/2) |]l0| , |l8| . 
The full line in Fig. 0(b), which plots our results for the 
same temperature for the insulating case, (n) = 1.0, in- 



2 



dicates the similarities between the 7% doped situation 
and the antiferromagnetic band structure at half-filling. 
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FIG. 2. The QMC spectral weight, of the 8 x 8 

Hubbard model with U = 8t as a function of w and fc: (a) at 
temperature T = 0.33t and 5% doping ((n) = 0.95), (b) at 
T = 0.25i and 7% doping ((n) = 0.93) and, (c) at T = 0.33t 
and 28% doping. Dark (white) areas correspond to large 
(small) spectral weight. The full lines in (a) and (c) are 
tight-binding fits to the QMC data, the full line in (b) denotes 
the U — 8t QMC result for (n) = 1.0, i.e. the insulating case. 

The physical origin of a possible manifestation of in- 
sulating band features in the under-doped regime seems 
to be that around the temperature T* ~ 0.3i the antifer- 
romagnetic correlation length £ becomes larger than the 
lattice spacing clq, i.e. for T — 0.25* we find £ ~ 1.2 > 1, 
whereas for T = 0.33i we find £ ~ 0.5 < 1 @. The de- 
cisive role played by antiferromagnetic spin fluctuations 
for the pseudogap and the under-doped band structure 
is substantiated by the behavior of Xs c(q, w), shown in 
Fig. 1 
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FIG. 3. Two-particle spin (uj sp in) and charge (uJ c harge 
in the inset) excitations as function of k, (a) under-doped 
((n) = 0.91), (b) over-doped ((n) = 0.72). Note the drastic 
change in the spin response as function of doping. 

Consider first Fig. ||(a), which exhibits the spin ex- 
citation energy as function of momentum for the 8x8 
Hubbard model with U = 8t and T = 0.33i at a doping of 
9% away from half-filling. Here, dark (white) areas again 
correspond to a large (small) spectral weight for spin ex- 
citations and the full line gives the spin- wave dispersion, 
i.e. u SDW (q) = 2J(1 - {ef/tt) 2 ) 1 / 2 , where e$ denotes 
the tight-binding energy e$ = —2t(cos(q x ) + cos(q y )). 
This figure reveals that in the under-doped regime up 
to maximal doping (13%) the spin excitations have an 
energy dispersion that still closely follows the antiferro- 
magnetic spin-wave dispersion. The main weight is con- 
fined to Q — (jr, 7r) with a small spread AQ and to a 
small but finite energy oj(Q) = to*. This sets an energy 
scale J for spin excitations. This result is in accordance 
with recent neutron-scattering data |2~fj(| and also with 
other experiments that have directly shown that prop- 
agating spin-waves with energies of order J still exist 
even at optimal doping pi] . When the temperature is 
lowered below T* ~ 0.3t, AQ slightly decreases and the 
antiferromagnetic correlation length £(T) (~ 1/AQ) be- 
comes larger than the lattice spacing. As a consequence, 
the single-particle hopping is now strongly renormalizcd 
by the short-range antiferromagnetic order resulting in a 
bandwidth (Fig. ||(b)) also of order of (a few) J. This 
renormalization is strongest at (tt, 0) and, thus, directly 
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responsible for the pseudogap: because of the "match- 
ing" condition, (tt,0) ± Q = (0,7r), the bare photohole 
created by ARPES for k — (n, 0) couples strongly to two- 
particle excitations whose spectral function peaks near 
momentum Q = (tt,it) 0. Using a standard diagram- 
matic evaluation of the self-energy like the FLEX sum- 
mation p2[ , and approximating the peak in Im Xs (<f, w) 
by Im Xs ~ S(q — Q)[—S(uj — uj*) + 6(u) +ui*)], one obtains 
for id < —id* < the self-consistency equation 



A(k,uj) 



U 2 A(k-Q,LU+LU* 



2 2 

u-E?-ReS(A,w)) + Im£(fc,w) 



, (2) 



since Im£(fc, ui) ~ —U 2 A(k — Q,ui + to*). As a conse- 
quence of this self-consistency requirement, multiple, re- 
peated spin- wave excitations ( "shake ups" ) accompany 
the bare carrier motion, lead to incoherent contributions 
to the electronic spectrum at (jr, 0) and spread spec- 
tral weight to lower energies. In the over-doped regime, 
shown in Fig. [|(b), the sharpness both in q and to in 
Im x s (g, w), and thus the "matching" condition is lost 
completely and the coupling to the bare photohole is 
weak and broad in energy. The crossover to a new energy 
scale (t rather than J) is also reflected in a corresponding 
crossover in the single-particle bandwidth (~ 8t) in Fig. 
§(c). In contrast to the spin excitations, the charge exci- 
tations are already broad and structureless in the under- 
doped regime, with the energy spreading essentially over 
the non-interacting bandwidth ~ 8t (inset in Fig. ||(a)) 
II- 

In this letter, we have addressed the microscopic mech- 
anism behind the opening of a pseudogap and its evolu- 
tion from under-doped to over-doped regimes. Our key 
results are (i) that the dynamical QMC results for the 
two-dimensional Hubbard model reproduce both the mo- 
mentum (i.e. the d x 2 _ y 2 -symmetry) and doping depen- 
dence of the pseudogap and (ii) that it is the unique 
doping dependence of the magnetic response which es- 
tablishes the first numerical proof that the pseudogap is 
due to antiferromagnetic spin correlations in the under- 
doped regime. 
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